clear


**Table 1
use "data_for_summary_statistics.dta", clear

foreach var of varlist head_male_hh head_age_hh head_ed_hh mother_age respondent_education fasting_mom age_child_months hhsize_hh TLU goodroof electricity farm_size poorhouse_hh distance_20k log_distance_20k  dds_mom_2019 dds_mom_2nd mindds_mom dds_child_2019  food_gap food_insecure staples_beans_nuts animal_source veg_fruits region1 region2 region3 region4 interview_weekend interview_time interview_afternoon call_attempt_num{
	capture ttest `var' if round==3, by(treat_all_round) level(95)
	gen `var'_mco = r(mu_1)
	gen `var'_mtr = r(mu_2)
	gen `var'_mdf = r(mu_1)-r(mu_2)
	gen `var'_li = (r(mu_1)-r(mu_2)) + r(se)*invt(r(df_t), 0.025)
	gen `var'_ui = (r(mu_1)-r(mu_2)) + r(se)*invt(r(df_t), 0.975)
}
keep head_male_hh_mco - call_attempt_num_ui
collapse (first) head_male_hh_mco - call_attempt_num_ui

gen str6 vars="mean", before(head_male_hh_mco)

preserve 
keep vars *_li
renvars head_male_hh_li head_age_hh_li head_ed_hh_li mother_age_li respondent_education_li fasting_mom_li age_child_months_li hhsize_hh_li TLU_li goodroof_li electricity_li farm_size_li poorhouse_hh_li distance_20k_li log_distance_20k_li dds_mom_2019_li dds_mom_2nd_li mindds_mom_li dds_child_2019_li food_gap_li food_insecure_li staples_beans_nuts_li animal_source_li veg_fruits_li region1_li region2_li region3_li region4_li interview_weekend_li interview_time_li interview_afternoon_li call_attempt_num_li \ head_male_hh head_age_hh head_ed_hh mother_age respondent_education fasting_mom age_child_months hhsize_hh TLU goodroof electricity farm_size poorhouse_hh distance_20k log_distance_20k  dds_mom_2019 dds_mom_2nd mindds_mom dds_child_2019  food_gap food_insecure staples_beans_nuts animal_source veg_fruits region1 region2 region3 region4 interview_weekend interview_time interview_afternoon call_attempt_num
tempfile lower_bound
save `lower_bound'
restore

preserve
keep vars *_ui
replace vars="upper" 
renvars head_male_hh_ui head_age_hh_ui head_ed_hh_ui mother_age_ui respondent_education_ui fasting_mom_ui age_child_months_ui hhsize_hh_ui TLU_ui goodroof_ui electricity_ui farm_size_ui poorhouse_hh_ui distance_20k_ui log_distance_20k_ui dds_mom_2019_ui dds_mom_2nd_ui mindds_mom_ui dds_child_2019_ui food_gap_ui food_insecure_ui staples_beans_nuts_ui animal_source_ui veg_fruits_ui region1_ui region2_ui region3_ui region4_ui interview_weekend_ui interview_time_ui interview_afternoon_ui call_attempt_num_ui \ head_male_hh head_age_hh head_ed_hh mother_age respondent_education fasting_mom age_child_months hhsize_hh TLU goodroof electricity farm_size poorhouse_hh distance_20k log_distance_20k  dds_mom_2019 dds_mom_2nd mindds_mom dds_child_2019  food_gap food_insecure staples_beans_nuts animal_source veg_fruits region1 region2 region3 region4 interview_weekend interview_time interview_afternoon call_attempt_num
tempfile upper_bound
save `upper_bound'
restore 

preserve
keep vars *_mtr
replace vars="mn_treat" 
renvars head_male_hh_mtr head_age_hh_mtr head_ed_hh_mtr mother_age_mtr respondent_education_mtr fasting_mom_mtr age_child_months_mtr hhsize_hh_mtr TLU_mtr goodroof_mtr electricity_mtr farm_size_mtr poorhouse_hh_mtr distance_20k_mtr log_distance_20k_mtr dds_mom_2019_mtr dds_mom_2nd_mtr mindds_mom_mtr dds_child_2019_mtr food_gap_mtr food_insecure_mtr staples_beans_nuts_mtr animal_source_mtr veg_fruits_mtr region1_mtr region2_mtr region3_mtr region4_mtr interview_weekend_mtr interview_time_mtr interview_afternoon_mtr call_attempt_num_mtr \ head_male_hh head_age_hh head_ed_hh mother_age respondent_education fasting_mom age_child_months hhsize_hh TLU goodroof electricity farm_size poorhouse_hh distance_20k log_distance_20k  dds_mom_2019 dds_mom_2nd mindds_mom dds_child_2019  food_gap food_insecure staples_beans_nuts animal_source veg_fruits region1 region2 region3 region4 interview_weekend interview_time interview_afternoon call_attempt_num
tempfile mn_treat
save `mn_treat'
restore 

preserve
keep vars *_mdf
replace vars="mn_diff" 
renvars head_male_hh_mdf head_age_hh_mdf head_ed_hh_mdf mother_age_mdf respondent_education_mdf fasting_mom_mdf age_child_months_mdf hhsize_hh_mdf TLU_mdf goodroof_mdf electricity_mdf farm_size_mdf poorhouse_hh_mdf distance_20k_mdf log_distance_20k_mdf dds_mom_2019_mdf dds_mom_2nd_mdf mindds_mom_mdf dds_child_2019_mdf food_gap_mdf food_insecure_mdf staples_beans_nuts_mdf animal_source_mdf veg_fruits_mdf region1_mdf region2_mdf region3_mdf region4_mdf interview_weekend_mdf interview_time_mdf interview_afternoon_mdf call_attempt_num_mdf \ head_male_hh head_age_hh head_ed_hh mother_age respondent_education fasting_mom age_child_months hhsize_hh TLU goodroof electricity farm_size poorhouse_hh distance_20k log_distance_20k  dds_mom_2019 dds_mom_2nd mindds_mom dds_child_2019  food_gap food_insecure staples_beans_nuts animal_source veg_fruits region1 region2 region3 region4 interview_weekend interview_time interview_afternoon call_attempt_num
tempfile mn_diff
save `mn_diff'
restore 

keep vars *_mco
renvars head_male_hh_mco head_age_hh_mco head_ed_hh_mco mother_age_mco respondent_education_mco fasting_mom_mco age_child_months_mco hhsize_hh_mco TLU_mco goodroof_mco electricity_mco farm_size_mco poorhouse_hh_mco distance_20k_mco log_distance_20k_mco dds_mom_2019_mco dds_mom_2nd_mco mindds_mom_mco dds_child_2019_mco food_gap_mco food_insecure_mco staples_beans_nuts_mco animal_source_mco veg_fruits_mco region1_mco region2_mco region3_mco region4_mco interview_weekend_mco interview_time_mco interview_afternoon_mco call_attempt_num_mco \ head_male_hh head_age_hh head_ed_hh mother_age respondent_education fasting_mom age_child_months hhsize_hh TLU goodroof electricity farm_size poorhouse_hh distance_20k log_distance_20k  dds_mom_2019 dds_mom_2nd mindds_mom dds_child_2019  food_gap food_insecure staples_beans_nuts animal_source veg_fruits region1 region2 region3 region4 interview_weekend interview_time interview_afternoon call_attempt_num 

append using `mn_treat'
append using `mn_diff'
append using `lower_bound'
append using `upper_bound'

xpose, clear

rename v1 mean_control
rename v2 mean_treatment
rename v3 mean_difference
rename v4 lower_bound
rename v5 upper_bound

drop in 1
gen str6 var_list=".", before(mean_control)

replace var_list="head_male_hh" in 1
replace var_list="head_age_hh" in 2 
replace var_list="head_ed_hh" in 3 
replace var_list="mother_age" in 4 
replace var_list="respondent_education" in 5
replace var_list="fasting_mom" in 6
replace var_list="age_child_months" in 7 
replace var_list="hhsize_hh" in 8 
replace var_list="TLU" in 9 
replace var_list="goodroof" in 10 
replace var_list="electricity" in 11 
replace var_list="farm_size" in 12 
replace var_list="poorhouse_hh" in 13 
replace var_list="distance_20k" in 14 
replace var_list="log_distance_20k" in 15  
replace var_list="dds_mom_2019" in 16 
replace var_list="dds_mom_2nd" in 17 
replace var_list="mindds_mom" in 18  
replace var_list="dds_child_2019" in 19  
replace var_list="food_gap" in 20  
replace var_list="food_insecure" in 21  
replace var_list="staples_beans_nuts" in 22  
replace var_list="animal_source" in 23  
replace var_list="veg_fruits" in 24
replace var_list="region1" in 25 
replace var_list="region2" in 26 
replace var_list="region3" in 27 
replace var_list="region4" in 28 
replace var_list="interview_weekend" in 29 
replace var_list="interview_time" in 30 
replace var_list="interview_afternoon" in 31 
replace var_list="call_attempt_num" in 32

format mean_control mean_treatment mean_difference lower_bound upper_bound %9.2f



**Supplementary table 1
use "data_summary_statistics.dta", clear

reg treat_all_round head_male_hh head_age_hh head_ed_hh mother_age respondent_education fasting_mom age_child_months  hhsize_hh TLU goodroof electricity farm_size poorhouse_hh distance_20k log_distance_20k dds_mom_2019 dds_mom_2nd dds_child_2019 food_gap food_insecure region1 region2 region3 region4 interview_weekend interview_time interview_afternoon call_attempt_num if round==2

estimate store m1

estout m1, cells(b(star fmt(3)) ci(par fmt(2))) starlevels( * 0.10 ** 0.05 *** 0.010) stat(F r2 N, fmt(2 2 0))
esttab m1 using "Join_significant_test.rtf", b(%9.3f) ci(%9.2f) star(* 0.10 ** 0.05 *** 0.01) stat(F r2 N, fmt(2 2 0)) coeflabels(_cons constant) title("Supplementary Table 1: Joint significance test") replace 
